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Abstract 

It is shown that Uhlmann's parallel transport of purifications along a path 
of mixed states represented by 2 x 2 density matrices is just the path ordered 
product of Thomas rotations. These rotations are invariant under hyperbolic 
translations inside the Bloch sphere that can be regarded as the Poincare ball 
model of hyperbolic geometry. A general expression for the mixed state geometric 
phase for an arbitrary geodesic triangle in terms of the Bures fidelities is derived. 
The formula gives back the solid angle result well-known from studies of the pure 
state geometric phase. It is also shown that this mixed state anholonomy can be 
reinterpreted as the pure state non-Abelian anholonomy of entangled states living 
in a suitable restriction of the quaternionic Hopf bundle. In this picture Uhlmann's 
parallel transport is just Pancharatnam transport of quaternionic spinors. 

I. Introduction 



Mixed state geometric phases as introduced by Uhlmann [1] provide a natural 
generalization of the well-known geometric phases [2] characterizing the geometric 
properties of unitarily or nonunitarily evolving pure states. Recently this idea 
of mixed state anholonomy was reconsidered within an interferometric approach 
[3,4], along with an alternative formulation of mixed state phases [5] . Possible 
experiments for confirming the appearance of such phases has been proposed and 
conducted (see e.g. [6]) and their relevance in the evolution of systems subjected 
to decoherence through a quantum jump approach has been stressed [7]. 



The simplest example of an evolving one-qubit system giving rise to a path 
in the space of nondegenerate 2x2 density matrices was studied by many authors. 
Uhlmann himself established a formula for the geometric phase for geodesic tri- 
angles and quadrangles drawn on the surface of a sphere of constant radius inside 
the Bloch ball B [8]. An explicit formula for the anholonomy along an arbitrary 
geodesic segment in B with respect to the Bures metric was presented by Hiibner 
[9] . In a recent paper for evolving systems giving rise to geodesic triangles defined 
by Bloch vectors in B of fixed magnitude, Slater [10] compared Uhlmann's geomet- 
ric phase with the interferometric approach of [5] . The aim of the present paper is 
to point out for such systems an interesting connection between Uhlmann's parallel 
transport and the phenomenon of Thomas precession. Using this correspondence 
we present a formula valid for an arbitrary geodesic triangle inside B. Our method 
is motivated by previous observations of Ungar [11] that hyperbolic geometry can 
be useful in describing the physical and mathematical phenomena associated with 
one qubit density matrices. 

The organization of the paper is as follows. In Section II. we briefly recall 
the background material needed for the definition of Uhlmann's anholonomy for 
mixed states. In Section III. using the hyperbolic parametrization of one qubit 
density matrices we show that Uhlmann's parallel transport can be expressed as 
the path ordered product of suitably defined Thomas rotations arising from the 
multiplication of two hyperbolic rotations (Lorentz-boosts) . Here the invariance of 
these rotations under the so called hyperbolic translates in the Poincare ball model 
of hyperbolic geometry is also established. In Section IV. we present our anholon- 
omy formula valid for an arbitrary geodesic triangle in the space of nondegenerate 
one qubit density matrices. Our result for the pure state limit gives back the solid 
angle rule well-known from studies concerning the geometric phase [2,5,10]. Here 
an error of Ref. [10] is also pointed out. Section V. is devoted to establishing a 
connection between our mixed state anholonomy and the pure state non-Abelian 
geometric phase. Section VI. is left for the comments and conclusions. 
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II. Mixed state anholonomy 



According to Uhlmann [1] mixed state anholonomy can be denned by lifting 
the curve pit) living in the space of strictly positive density operators to the space 
of its purifications such that the representative curve Wit) of these purifications is 
parallel translated with respect to a suitable connection. W purifies p if we have 
p = WW^ and TrWW^ = 1. For the special case of purifications W that are 
elements of 7i <8> 7i* with 7i is a finite n-dimensional Hilbert space and give rise 
to strictly positive density matrices we have W G GL{n,C). It is obvious that 
the process of purification is ambiguous, W and WU with U G U(n) gives rise to 
the same p. Hence we have a (trivial) principal bundle with total space GL(n, C) 
base X>+ = GL(n, C)/U(n) and fiber U(n). According to the connection defined 
by Uhlmann two purifications W\ and W 2 giving rise to p\ and p 2 respectively are 
parallel iff 

w\w 2 = w\w x > 0. (1) 

1/2 1/2 

Using the polar decompositions Wi = p/ Ui and, W 2 = p 2 U 2 , that are arising 
as the right translates by U(n) of the global section p 1 / 2 one can check that U\ 
and U 2 related as 

v — tt rrt - 1 / 2 -!/2/ 1/2 1/2x1/2 / Q \ 

Y 21 = U 2 U{ = p 2 * Pl * (p/ p 2Pl > ) 1 (2) 

gives rise to purifications satisfying (1). Dividing our path p(t) G £>+ with < 
t < 1 into small segments one gets 



W 1 W^= \imX lsi X SlS2 ...X Sm0 p , X st = p; 1/2 {p] ,2 p s p] /2 ) 1 ' 2 p; 1/2 (3) 

where "lim" indicates the process of going to finer and finer subdivisions producing 
a continuous path C G Although Eq. (3) with the basic building blocks beeing 
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the X st was usually used in the literature, for later use we favour an alternative 
one expressed in terms of the Y st defined by (2) as 

W t Wi = lim p{ /2 Y lsi Y SlS2 . ..Y Sma pl /2 . (4) 

m^oo 

This expression (the anholonomy of the curve C) will be our basic one for the 
description of parallel transport of purifications over a path C G T>+ . For a closed 
path C we have po = pi hence for the trace of this expression we get 

Tr(W 1 Wt) = lim Tr(Y Ul Y SlS2 . . .Y Sm0 p ). (5) 

The quantity $ 3 = arg Tt(WiWq) is the generalization of the geometric phase for 
mixed states, and the magnitude v = |Tr(WiWg)| is the visibility [5]. 

III. Mixed state anholonomy as Thomas rotation 

Now we start discussing our main concern here, namely mixed state anholon- 
omy for a qubit system. For these systems p is an element of the interior of the 
usual Bloch-ball £>. Moreover, the space of purifications for strictly positive 2x2 
density matrices p G ~ Int£> is GL(2,C). Hence we have p = WW^ with 
W G GL(2, C). It is particularly instructive to regard this space as the space of 
normalized entangled states for a bipartite system i.e. to have H^H* ~ C 2 <g> C 2 , 
a description giving rise to our p upon taking the partial trace with respect to the 
second subsystem. 

For this we write the entangled state G C 2 <g> C 2 in the form 



m = J2 w ^\jk), ij^ijk®!/^, j,k = o,i w jk = ^=(^ c b d y (6) 

The normalization condition = 1 in this picture corresponds to the con- 

straint Tr(VFVF^) = 1. Moreover, calculating the trace of the pure state density 
matrix with respect to the second subsystem yields our p, i.e. in the 
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base we have p = Tr 2 |#)(#| = WW^ . Since we are considering strictly positive 
density matrices we have the constraint Detp 7^ 0. In terms of the matrix elements 
of W it means ad — be 7^ i.e. W G GL(2, C). It is well-known that for this bi- 
partite system the measure of entanglement is the concurrence C [12] which can 
be written as 



0<C=\ad-bc\<l. (7) 

Separable states corresponding to reduced density matrices with a zero eigenvalue 
are precisely the ones with C = 0. 

Let us parametrize our density matrix as p = |(/+ucr) = WW^ with |u| < 1. 
This means we have 2-u 3 = |a| 2 + |6| 2 — |c| 2 — |<i| 2 and u\-\-iu<i = ac+bd, and it is easy 
to check that C = ^/l — |u| 2 = 4Detp. The four real numbers —1 < ui, U2, us < 1 
and < U4 = C < 1 can be regarded as coordinates on the upper hemisphere of a 
three-dimensional sphere S 3 embedded in R 4 homeomorphic to Int£>. 

In the following it is convenient to introduce a new (hyperbolic) parametriza- 
tion for p by introducing the rapidities [13,14] 9 as 



|u| = tanh0 u , < 9 U < 00. (8) 

In this parametrization the concurrence is related to the quantity 7„ = cosh 9 U of 
special relativity as C = r )~ 1 . The reader can verify that in this case 



P 1/2 = \j\ (cosh|/ + sinh|u^y|L(^,u), u^^, (9) 

where L(9 U , u) is a Lorentz boost in the spinor representation. 

Using the relation M 2 - MTrM + DetM = (valid for 2 x 2 matrices) and 
its trace we have the formula 



M = , M + V5itg (10) 
V / TrM + 2 v / DetM 
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Now we calculate the quantity Y uv of (2) with the density matrices p u and p v . For 

1 /2 1 /2 

this we insert M = p</ /9 u pt/ in Eq. (10). First we note that the square of the 
denominator of this equation has the form 

Tr(p uPv ) + 2y/Det( PuPv ) = ^(1 + uv + C U C V ) = F(u, v), (11) 
where F(u, v) is the Bures fidelity. Moreover, by virtue of Eq. (9) its enumerator 

2/2 1/2 

multiplied by p u p v has the form 

pI/ 2 pI /2 + ^pu 1/2 Pv 1/2 = \ Vc^c~v (L(e u , u)L(e v , v) + L(e u , -u)l(9 v , -v)) . 

(12) 

Recalling that the composition of two boosts can be written as another boost times 
a Thomas rotation the right hand side of Eq. (12) can be written as 

Q 

L(6 W , w)R(a, n) + L(6 W , -w)R(a, n) = 2 cosh -^R(a, n) (13) 

where R{a, n) = cos ^1 + i sin fan is the Thomas rotation matrix in the spinor 
representation. By writing out explicitly the product of two boosts in the spinor 
representation the dependence of the quantities a , ^, w,non the original ones 
9 U , 6 V and u, v can be established (see e.g. [15] and references therein). One of such 
formulas we need is the hyperbolic law of cosines [13,15] which can be written as 

cosh 9 W = cosh 9 U cosh 9 V + sinh 9 U sinh 6*„uv = (1 + uv) . (14) 

Using this we have 2 cosh 6 -f = 2^^^-. Putting this into Eqs. (12-13) and 
using Eq. (11) we obtain our result 

Y uv = RZ^PZ^ipV 2 PuPl /2 ) 1/2 = R{«, n). (15) 
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Hence according to Eq. (4) Uhlmann's parallel transport can be understood as a 
sequence of Thomas rotations. 

For the calculation of the explicit form of R(a, n) we use the left hand side of 

1/2 1/2 

Eq. (12) and the explicit forms of p u and p v obtained from Eq. (9) by expressing 
the hyperbolic functions in terms of C u , and C v . 

A/ ' 2= 27TTK {1 + c " + ua) ' py2 =viW 1 + c " +w) - (16) 

Collecting everything we get 



Y uv — R(a, n = cos —I + % sin — na = , = • 

2 2 ^[(l + CXl + GO+uvP + lux v| 2 

(17) 

(Compare this explicit formula with the implicit one of Eq. (2) of Ref. [9].) From 
this 



a u x v u x v o . 

tan— = - — — — , n=- r . 18 

2 (1 + C u )(l + C v ) + uv J |uxv| v ' 

According to [1] this parallel transport is the one occurring along the shortest 
geodesic with respect to the Bures metric between the two points u and v in the 
interior of the Bloch-ball. Since every smooth curve can be approximated by a 
sequence of geodesic segments Eq. (4) can be regarded as such an approximation. 
According to our result the parallel transport of purifications along a smooth curve 
in B can be represented as the path ordered pruduct of Thomas rotations. 

Let us examine the (17) expression for the anholonomy transformation ma- 
trix. First we introduce a new parametrization 



u 9 U A V 6y ^ . . 

tanh — u, b = — = tanh — v. 19 



l + C u 2 1 + C V 2 
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It is clear that a and b are still elements of the Bloch-ball, they are of the same 
direction but different length. In terms of these new variables R(a, n) can be 
written as 



R(a, n) = R(a, b) = 1 J =L. G Stf(2), (20) 

a/1 + 2ab + |a| 2 |b| 2 

where we have used the formula |a x b| 2 = |a| 2 |b| 2 — (ab) 2 . Notice that in this 
notation the Bures fidelity is related to the square of the denominator of this 
formula via the identity 



F ( u > v) = i(l +C U )(1 +C)(1 + 2ab+ |a| 2 |b| 2 ) = 1 - {1 + l^ l +lhl2 y (21) 

Here the second equality of Eq.(21) also reveals the relationship of < F(u, v) < 1 
to the distance on the upper hemisphere of S 3 , as can be checked by stereographic 
projection from the south pole of S 3 to R 3 that maps the upper hemisphere of S 3 
to B. An alternative form of (21) is F(u, v) = cos 2 y, where < A < n is the 
geodesic distance between a and b with respect to the metric on B arising via this 
stereographic projection. 

Now let us define the hyperbolic-translation [16] of the vector b by the vector 

a as 



r a (b) = 



(1-1 


a 


2 )b + (1 + 2ab+ |b| 2 )a 


1 + 2ab + 


a| 2 |b; 
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(22) 



Denoting a' = r a (b) one can show that |a'| 2 = |a + b| 2 /(l + 2ab+ |a| 2 |b| 2 ) < 1 i.e. 
this transformation maps B onto itself. Using this we get (1 + 2a'b + |a'| 2 |b| 2 ) (1 + 
2ab + |a| 2 |b| 2 ) = (1 + 2ab + |b| 2 ) 2 . Now it is easy to establish the formula 



i?(a, b) = fl(r a (b), b) = R{a, r b (a)). (23) 

Equation (23) states the important result that the anholonomy for Uhlmann's par- 
allel transport is invariant with respect to hyperbolic-translations of the Bloch-ball 
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regarded as the Poincare ball modell of hyperbolic geometry. These properties were 
called by Ungar [11] "left-loop" and "right- loop" properties in his study of density 
matrices and gyrovector spaces. In this way we established an implementation 
of his abstract setting up on Uhlmann's parallel transport of purifications. We 
note in closing that the (23) properties can be used to find deformations of curves 
consisting of geodesic segments having the same anholonomy, a property that can 
be useful for the experimental verification of Uhlmann parallelism in this most 
general setting up. 

IV. Geodesic triangles 

In this section we use the results of the previous section to derive an explicit 
formula for the anholonomy of a special closed path: the geodesic triangle. Note, 
that this problem has already been considered in Refs. [8] and [10] for the three 
points of the triangle lying on a spherical shell of B of fixed concurrence. Here 
we consider the general case and chose three arbitrary points u, v and w in the 
interior of the Bloch-ball. The corresponding concurrences are C u , C v and C w . 
We renormalize our vectors u, v and w as in Eq. (19), the resulting vectors still 
belonging to Int£> are a, b and c. According to Eqs. (4,17,20) in order to calculate 
the anholonomy we have to evaluate the quantity 

5 S 

7£(a, b, c) = R(a, c)R(c, b)R(b, a) = cos - + % sin -ma, (24) 

where the last equality expresses the fact that the resulting matrix should also 
have an SU(2) form. 

By virtue of Eqs. (20-21) we can extract a factor from 1Z(&, b, c) of the form 
|(1 + C u )(l + C v )(l + C w )/y/F(u, v)F(w, v)F(v, u). Hence we merely have to 
evaluate the matrix 



J = (I+a-c)(I+c-b)(I+b-a), where a-c = (acr)(ca) = (ac)/+z(ax c)a etc. 

(25) 
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Straightforward calculation yields the result 



J = (1 + aW + (1 + a 2 )(bc) + (1 + 6 2 )(ca) + (1 + c 2 )(ab)) I 



+ [c- (ca)a, b - (ab)a] + \ ((1 - a 2 ) [6, c] - (1 - 6 2 ) [c, a] - (1 - c 2 ) [a, 6]) , (26) 



where a = (acr), a 2 = |a| 2 etc. and [a, c] denotes the commutator of the corre- 
sponding matrices. By virtue of the relations a 2 = (1 — C u ) /(l + C u ) etc. and the 
(11) definition of the Bures fidelity we obtain for that part of 7£(a, b, c) which is 
proportional to the identity matrix the formula 



cos 6 = F(u,w) + F(w,v) + F(v,u)-1 
2 2 v /F(u,w)F(w,v)F(v,u) 

In order to also find the axis of rotation we introduce the vectors 



_ (l + a 2 )c- (l + 2ac-c 2 )a _ (1 + a 2 )b - (1 + 2ab - 6 2 )a 
P= 1 + 2ac + a 2 c 2 ' Q= 1 + 2ab + a 2 b 2 ' ^ ' 

With these vectors it is straightforward to check that 



^(a,b,c)= P + P^ + '(P^ , (29) 
a/1 + 2pq + p z q^ 

hence the angle 5 and axis m of the resulting Thomas rotation is given by the 
expressions 

<5 |p x q| P x Q 

tan- = t^- — — ^t-, m=-^- — 30 

2 (1 + pq) |pxq| 
A comparison of Eq. (28) and (22) shows that the formulae for p and q 
up to some crucial sign changes look like the hyperbolic translates, p and q are 
"translates" by a of c and b. However, by virtue of the relations p 2 = |a — c| 2 / (1 + 
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2ac + a 2 c 2 ), and q 2 = |a — b| 2 / (1 + 2ab + a 2 b 2 ) one can see that these "spherical 
translates" are not mapping B (homeomorphic to the upper half of S 3 ) onto itself. 
These are rather isometries of the full S 3 with its metric given by Eq. (21). 

In order to gain some more insight into the geometric meaning of the (28) 
" spherical translate" we notice that 



|pp |c — a 

Now recall [16] the definition of the transformation cr^ : R 3 — > R 3 

^W = a+-^(x-a), (32) 

which is the inversion with respect to a sphere S 2 in R 3 centered at a with radius 
r. It is now obvious that (31) is an inversion of c with respect to a sphere centered 
at a with radius r 2 = 1 + a 2 . Moreover, since the transformation p — > p/p 2 is also 
an inversion with respect to the sphere centered at the origin with radius 1, we 
obtain the following result, p (resp. q) is the result of two inversions applied to 
the point c (resp. b). One of the inversions is defined by the point a, the point 
we have chosen as the starting one for the traversal of the geodesic triangle. 

Now let us calculate the geometric phase corresponding to our geodesic tri- 
angle! First we notice that 



9 9/ 9^9 l + 2bc + 6 2 C 2 

l + 2pq + pV = 1 + a 2 2 - — ^-r— 5-^. 33 

v ; (l + 2ab + a 2 6 2 )(l + 2ac + a 2 c 2 ) 

From Eqs. (5) and (24) we get Tr(7?.(u, v, w)p u ) = cos | + i sin |(um) hence after 
recalling Eqs. (19), (21), (27-30) and (33) straightforward calculation yields the 
result 



;a , F (u,v ) + F v,w ) + F w,u - 1 - k v x w , 

ve i* s = Tr(n Pu ) = K ' V ' V ' 2 -. 34 

2 v /F(u,v)F(v,w)F(w,u) 

Hence our formula for the mixed state geometric phase takes the form 
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ill(v X w) 

tan& a = ——, r —, r —, r . (35) 

9 F(u,v) + F(v,w) + F(w,u) - 1 v 1 

In the pure state limit we have C u = C v = C w = 0. Since the vectors u, v and w 
are of unit norm we denote them in this special case as ni, n 2 and n3. Recalling 
Eq. (11) for the Bures- fidelity from (35) we obtain the result 



+ * + ni(n 2 x n 3 ) 

tan<P„ = — tan— = (36) 

2 1 + nin 2 + n 2 n 3 + n 3 ni 

which is the formula for the tangent of minus half the solid angle [17] corresponding 



to the geodesic triangle on the surface of the unit sphere S 2 i.e. we have $ g = 



£2 
2 • 



Note that a more familiar form for Q is given by [17] 



O cos 2 %i + cos 2 ^f- + cos 2 %i - 1 1 + nin 2 + n 2 n 3 + n 3 ni 

2 2 co^co^3co^i V2(l + n 1 n 2 )(l + n 2 n 3 )(l + n 3 n 1 )' 

(37) 

with n^iij = cos9ij. Comparing this with the pure state limit of Eq. (27) we 
see that in this case the Thomas rotation angle is just the solid angle, i.e. 5 = 0. 
Hence for pure states we get back to the results well-known from studies concerning 
the ordinary geometric phase. 

Eq. (34) is the most general formula that defines the visibility v and the 
Uhlmann mixed state geometric phase $ 5 valid for an arbitrary geodesic triangle 
defined by the points u, v, w inside the Bloch ball B. The usual geometric phase 
is obtained in the limiting case of sending all of the points to the boundary of 
B representing pure states. As a further investigation of formula (34) let us now 
consider the important special case studied by Slater [10] when |u| = |v| = |w| = 
r = const! Let u = rni, v = rn 2 and w = rn 3 , then F(u, v) = 1 + ir 2 (nin 2 — 1) 
etc. By virtue of (35) we have the result 



r 3 ni(n 2 x n 3 ) r 3 /i O 

— tctn 

4(1 - r 2 ) + r 2 (l + nin 2 + n 2 n 3 + n 3 ni) 4(l-r 2 )+r 2 ^ 2' 

(38) 
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where \i = 1 + 11^2 + n 2 n 3 + n 3 n! a notation used in Ref [10]. This result is in 
contrast vith the claim of Slater (see Eq. (18) of Ref. [10]) 



tan$ Slater = _ Q 

9 4 + (p - 10)r 2 + 6r 4 2 V ; 

Notice that after the replacement 6r 4 — > 6r 2 his result would reproduce the correct 
one of Eq.(38). 

Using ideas of interferometry an alternative definition (different from the one 
as given by Uhlmann) for the mixed state geometric phase appeared in Ref. [5]. 
In this approach the result for the situation studied above is [5] 

tan^ nt = -rtan— . (40) 

From Eqs. (38) and (40) we see that the ratio tan $ 5 / tan $^ nt is r 2 / u/r 2 / u+4(l-r 2 ) 
hence the the two different types of mixed state phases are equal merely in the 
pure state (r = 1) case. (In Ref. [10] it was claimed that the two phases are equal 
also for the nontrivial case with r = a/2/3, a possibility clearly following from 
the erroneous result of Eq. (39).) The fact that the two approaches give different 
results for the mixed state anholonomy was first stressed in Refs. [10] and [4]. 

Closing this section we check that the formula for the visibility (i.e. the 
magnitude of the rhs. of Eq. (34) ) gives the result v = 1 in the pure state limit. 
The visibility is 



(F(u, v) + F(v, w) + F(w, u) - l) 2 + lU 2 
' ~~ V 4F(u,v)F(v,w)F(w,u) ' 1 ' 

where V = u(v x w) is the volume of the parallelepiped spanned by the triple u, 
v and w . In the pure state limit we get 



" 2+y2 (42) 



2(1 + nin 2 )(l + n 2 n 3 )(l + n 3 ni) 
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Using the First Law of Cosines [16] cos 6>i 2 — cos #23 cos #31 = sin #23 sin 6 31 cos 7 
of spherical trigonometry, where 7 is the angle of the spherical triangle at the 
point defined by n 3 , and the relations rijiij = cos 0y and sin 6*23 sin 6*31 sin 7 = V 
straightforward calculation yields the expected result v = 1. Hence, for the mixed 
state case 5 is different from Q and v 7^ 1 (um 7^ 1) properties also shown by 
analysing the alternative formula v = \J cos 2 8/2 + sin 2 5/2 (um) 2 to be compared 
for r = constant with formula (26) of Ref [5]. 

V. Mixed state anholonomy and quaternionic phases 

Representing the space of purifications as the space of normalized entangled 
states in C 2 <g> C 2 (see Eq. (6)) we have the possibility to build up a dictionary 
between the nomenclatures of the mixed and the pure state anholonomies. In 
order to do this recall that due to the constraint (^l^) = TrWWt = 1 the space 
of such purifications is the seven sphere S 7 . Let us parametrize the matrix W of 
(6) as 

W = -^=(Qo + iQi), Qo = a I - iajo-j, Qi = (3 Q I - ifijcrj. (43) 

(Summation over repeated indices is understood, <jj j = 1, 2, 3 are the Pauli matri- 
ces.) Notice that Eq. (43) amounts to a change of parametrization from the four 
complex numbers a, 6, c, d of Eq. (6) to the eight real ones and (3^ n = 0, 1, 2, 3. 
Explicitly we have 

a = (a + (3 3 ) + i(P - a 3 ), b = - a 2 ) - i{a x + (3 2 ) 

c= (p x + a 2 ) -i{on -#>), d= (a - /3 3 ) +i(/3 + a 3 ). (44) 

and 
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2«o = + did, 2a i = —36 — 3c, 2a 2 = 9?c — 3?6, 2«3 = 3d — 3a 

2/3 = 3a + 3d, 2/3i = m + Kc, 2/3 2 = 3c - 36, 2/3 3 = »a - 3«d, (45) 

where the symbols 3? and 3 refer to the real and imaginary parts of the corre- 
sponding complex numbers. Notice moreover, that the correspondence 



defines a mapping between a W G GL(2, C) and a quaternionic spinor (go, qi) T G 
H 2 i.e. we have the correspondence 



On the space of two component quaternionic spinors we can define an inner product 
(|) : H 2 x H 2 — > H as (q\p) = qoPo + QiPii i- e - we have quaternionic conjugation 
in the first factor. From the normalization condition TrVFVF^ = 1 it follows that 
+ PuPn = 1) i- e - the spinor (go,gi) T is normalized, ||g|| 2 = (q\q) = 1. It 
means that ct M and /3 M are Cartesian coordinates for the seven-sphere S 7 . 

Let us express our reduced density matrix p = Tt2\^){^\ = WW^ in terms 
of the matrices Qq and Qi corresponding to the quaternions qo and gj By virtue 
of (43) we have 



where we have used the normalization condition and the fact that the matrix 
QiQo — QoQi is an anti-Hermitian 2x2 one hence it can be expanded as —iajUj, 
j = 1, 2, 3 with tti, U2, us are real parameters of the one-qubit density matrix. For 



1 <-> — Z<Ti, j <-> — ZO"2, k <-> — ZO"3 



(46) 




go = a + aii + a 2 j + a 3 k, g x = /3 + /5 x i + (3 2 j + /9 3 k. (47) 




(48) 
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later use we also define the quantity uq = QiQ Q + QoQ\ which is two times the 
Hermitian part of QiQq- The Hermitian and anti-Hermitian parts of the matrix 
Q1Q0 correspond to the real and imaginary parts of the corresponding quaternion 
qiq hence we can define the quaternion 

u = uq + ui'i + u 2 j + u 3 k = 2qi%, Re(tt) = -{u + u) = u , Im(u) = -{u - u). 

(49) 

Let us define one more quantity 



u 4 = \q \ 



|2 



\qi\ 2 = QoQt-QiQl (50) 



Recall also from Section III. that that the coordinates jl = 0, 1, 2, 3, 4 are related 
to the complex numbers a, 6, c, d as 114 + iuo = ad — be, u\ + iu2 = ac + bd and 
2us = \a\ 2 + \b\ 2 — |c| 2 — \d\ 2 . Hence the concurrence is just C = \fu\-\- Uq. It 
is straightforward to check that uu + u\ = Ujxup, = 1, hence up, E S 4 i.e. it is 
an element of the four dimensional sphere. As a result of Eqs. (49-50) one can 
define a map 7r : S 7 — > S 4 . Notice that according to the explicit form of this 
map the transformation (right multiplication of the quaternionic spinor with a 
unit quaternion) 

^ere ss = 1 (51) 

leaves the coordinates up invariant. Since unit quaternions correspond to elements 
of SU{2) ~ S 3 the projection it defines a fibration (the second Hopf fibration 
[18]) of S 7 with base S 4 and fiber S 3 . Reinterpreting our quaternionic spinors 
as entangled states it is straightforward to show that this SU(2) fiber degree of 
freedom corresponds to the possibility of making local unitary transformations I® 
S : S e SU(2) in the second subsystem. This idea of representing entanglement via 
the twisting of a nontrivial fiber bundle was initiated in [19] and further developed 
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in Refs. [20] and [21]. Here we merely need one result from Ref. [20]: an element 
|g) G S 7 can be parametrized by points of S 4 minus the south pole (SP) as 

\q)=( qo )= - 1 ( 1 + U4 )s = \u)s, 88=1, up, G S 4 — {SP}. (52) 

Eq. (52) is a local section of our bundle. There are no global sections (i.e. ex- 
pressions like (52) nonsingular over all of S 4 ) which is just another way of saying 
that the Hopf bundle is nontrivial i.e. S 7 ^ S 4 x S 3 . Of course we can define 
alternative sections that are singular at different points, the (52) choice is dictated 
by convenience. 

Consider now three quaternionic spinors \q) = \u)s, \p) = \v)x and |r) = \w)y 
s,x,y E SU(2) representing entangled states \^) , |$) and |%)! Notice that the 
notation indicates that the corresponding quaternionic spinors are parametrized 
by the vectors u/x , and w/x which are elements of the open neighbourhood 
S 4 - {SP}. 

As a next step we consider the trivial subbundle £ of the Hopf bundle defined 
by the conditions uq = 0, and u± = C > 0. £ is a fiber bundle with an S 3 
fiber over the submanifold M. of the upper half hemisphere of 5 4 defined by the 
aforementioned constraints. It is easy to see that M. is topologically the upper 
half hemisphere of a three sphere defined by the coordinates C, u\,U2, u% and can 
be identified with the interior of the Bloch-ball of reduced density matrices Int£>. 
For more details on the structure of the bundle £ that has already been studied 
in the context of Uhlmann's connection see Ref. [22]. Let us suppose that our 
spinors \q), \p) and \r) define global sections of £ of the (52) form. This means 
that we set the parameter values uo = v$ = wo in expressions like Eq. (52) to 
zero, and the ones ^4,^4, W4 to C u , C v and C w . 

It is now straightforward to check that the unit quaternion (v\u) /\(v\u)\ is 
just Y vu of Eq. (17). Moreover employing the notation of Eq. (19) equation (24) 
can be written in the following form 
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v< \ W) (r\p) (p q) _ 

\(q\r) \ \ (r\p) \ \ {p\q)\ 

where it is now understood that the left hand side is also regarded as a unit 
quaternion. 

Notice now that Eq. (52) is just the quaternionic analogue of the polar 
decomposition. Indeed according to Eq. (16), the spinor \u) corresponds to the 

1 /2 

matrix p u ' , and the unit quaternion s to the SU (2) part of the U (2) matrix S of 

1 /2 

the polar decomposition W = pj S. Since £7(2) ~ £7(1) x SU(2) we only have to 
account for a complex phase, but this is fixed by our choice uq = when restricting 
to the subbundle £. (Notice that according to Eqs. (6-7) 144+^0 = Ce lK = 2DetW 
where tan« = uq/u^. Hence the uq 7^ case amounts to multiplying our (6) 
entangled state by a £7(1) phase.) 

Now let us write Eq. (4) for the geodesic triangle in the following form 



W 1 = AW = p 1 J 2 n(u, V ,w)p- 1 / 2 W , where p u = W W^ (54) 

1 /2 

Since the polar decomposition Wq = p u S corresponds to the (52) section of the 
bundle £, we can write 

W 1 = Apl/ 2 S = pl/ 2 SS*K(u, v, w)S. (55) 

Using the notation \q') for the quaternionic representative of W\ we see that the 
quaternionic version of Eq. (55) is \q') = \q)~sJZs. By virtue of Eq. (53) Uhlmann's 
parallel transport in £ in the quaternionic representation can be written as 

, m , s. (<l\r} (r\p) (g|g) (56) 
l?> " k> l<9|r>|Krb>||<p|g>r (56) 

It is clear that for an arbitrary geodesic polygon Eq. (56) has to be multiplied 
from the right by extra quaternionic phase factors corresponding to transitions to 
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the new points of the polygon. The geodesic rule obtained in this way is the non- 
Abelian analogue of the well-known one obtained for filtering measurements in the 
context of the usual geometric phase [23], [24]. In this picture each polygon V is 
decomposed into a sequence of geodesic triangles. Each triangle Aj gives rise to a 
Thomas rotation of the (24) form with angle 5j and axis nx, . Since V in general is 
not a planar curve the rotations corresponding to different triangles have different 
axes. As a result the total rotation angle is not the sum of the component rotations 
as was in the Abelian case corresponding to the ordinary geometric phase. In this 
more general case we have to combine rotations with different axis resulting in 
the appearance of a path ordered product. Going to finer and finer subdivisions 
any smooth closed curve C can be approximated by a suitable polygon T. The 
resulting quaternionic phase can be written as the path ordered exponent Ve $c A 
where the s-u (2) -valued gauge-field can be written as 



1 udu 

A = lm(u\du) = -Im u = ui'i + U2J + t^k, -u 4 = C u . (57) 

2 l + «4 

As was remarked in Refs. [20,22] A is just the pull-back of the restriction of the 
canonical (instanton) connection on the quaternionic Hopf bundle to the bundle 
S with respect to the section s = 1 (see Eq. (52)). 

VI. Conclusions 



In this paper we investigated Uhlmann's parallel transport as applied to a 
qubit system. In spite of beeing the simplest and hence best studied example this 
system still shows nice geometric properties have not fully been appreciated by 
the physics community. Our paper was intended to fill in this gap by explicitly 
working out these missing interesting details. First we have shown that the very 
special features of the qubit system enable one to reinterpret Uhlmann's parallel 
transport as a sequence of Thomas rotations. We have also shown some interesting 
connections with hyperbolic geometry. In particular we proved that the finite 
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Thomas rotations are invariant with respect to hyperbolic translates of the interior 
of the Bloch ball regarded as the Poincare model of hyperbolic geometry (see Eq. 
(23)). These observations should not come as a surprise since Uhlmann's parallel 
transport has its origin in the underlying Bures geometry [1] of the Bloch ball £>, 
that has already been related to the Poncare metric in hyperbolic geometry [13,20], 
moreover it is easy to see [20] that the Bures metric is conformally equivalent to 
the standard Poincare one. 

In section IV. we derived an explicit formula for the SU(2) anholonomy 
matrix in the case of a geodesic triangle (Eqs. (28-30)). From this an expression 
in terms of the Bures fidelities for the mixed state geometric phase and the visibility 
was derived (Eq. (35) and (41)). These general results were shown to give back 
in the pure state limit the corresponding ones known from studies concerning the 
usual geometric phase. As far as the author knows these formulae in their full 
generality have not appeared in the literature yet. The geometric significance of 
these expressions were elaborated, and an error that appeared in Ref. [10] was 
corrected. 

In Section V. we managed to reformulate our results concerning the mixed 
state anholonomy in terms of the pure state non-Abelian one. The idea was to 
reinterpret the space of purifications as the Hilbert space for an entangled two- 
qubit system. This trick enabled us to recast Uhlmann's parallel transport in yet 
another form i.e. in the one of a sequence of quaternionic filtering measurements 
(Eq. (56)). By going to finer and finer subdivisions we have recovered Uhlmann's 
parallelism as the Wilson loop over a gauge field which is a suitable restriction of 
the well-known instanton connection. 

The advantage of this quaternionic formalism is clear: Uhlmann's parallel 
transport for one qubit density matrices in this representation is just the quater- 
nionic analogue of the usual Pancharatnam transport extensively used in studies 
concerning the geometric phase [2, 25]. In this language two entangled states |\&) 
and |$) regarded as purifications for one-qubit density matrices are "in phase" iff 
their quaternionic representatives \q) and \p) satisfy the constraint: (p\q) is real 
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and positive. It is easy to check that this constraint is equivalent to the one as given 
by Eq. (1). Moreover, this rule provides a nice way of defining the difference of 
these entangled states in their local unitary transformations corresponding to the 
second subsystem. Indeed, consider and |<&) as above and define their relative 



<* 


*> 







their relative SU(2) " (quaternionic) phase" as jg-gj. Since U(2) ~ U(l) x SU(2) 
this convention defines a relative £7(2) "phase" for our entangled states. When 
the entangled states in question have the same reduced density matrices this U(2) 
transformation corresponds to the possibility of the observer in the second subsys- 
tem to rotate the shared state to |<&), via his freedom to employ local unitary 
transformations. In the general case using this definition we can compare the local 
unitary transformations (corresponding to the second subsystem) of two entangled 
states with different reduced density matrices. 

It is clear that these results imply many interesting applications. Apart from 
studying the generalization of our results for nonsingular n x n density matrices 
via the use of the anholonomy defined by Uhlmann's connection on the trivial bun- 
dle GL(n,C)/U(n), there is also the interesting possibility of studying quantum 
gates defined by anholonomy transformations over the stratification manifold of 
entangled qudit systems. Though some of these issues have already been partly 
discussed [20,22] we hope to report some new results in a subsequent publication. 
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